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ABSTRACT 

PSR B 1259-63 is a 48 ms radio pulsar in a highly eccentric 3.4 year orbit with a Be star SS 
2883. Unpulsed y-ray, X-ray and radio emission components are observed from the binary 
system. It is likely that the collision of the pulsar wind with the anisotropic wind of the Be 
star plays a crucial role in the generation of the observed non-thermal emission. The 2007 
periastron passage was observed in unprecedented details with Suzaku, Swift, XMM-Newton 
and Chandra missions. We present here the results of this campaign and compare them with 
previous observations. With these data we are able, for the first time, to study the details of 
the spectral evolution of the source over a 2 months period of the passage of the pulsar close 
to the Be star. New data confirm the pre-periastron spectral hardening, with the photon index 
reaching a value smaller than 1.5, observed during a local flux minimum. If the observed 
X-ray emission is due to the inverse Compton (IC) losses of the 10 MeV electrons, then 
such a hard spectrum can be a result of Coulomb losses, or can be related to the existence 
of the low-energy cut-off in the electron spectrum. Alternatively, if the X-ray emission is 
a synchrotron emission of very high energy electrons, the observed hard spectrum can be 
explained if the high energy electrons are cooled by IC emission in Klein-Nishina regime. 
Unfortunately the lack of simultaneous data in the TeV energy band prevents us from making 
a definite conclusion on the nature of the observed spectral hardening and, therefore, on the 
origin of the X-ray emission. 

Key words: pulsars : individual: PSR B1259-63 - X-rays: binaries - X-rays: individual: 
PSR B 1259-63 



1 INTRODUCTION 

PSR B 1259-63 is a ~48 ms radio pulsar located in an eccentric 
(e~0 87), 3.4 year orbit with a Be star SS 2883 jjohnston et"aH 
Il992h . This system is kn own to be hig h ly var iable on an or- 
bital time scal e in radio (jjohnston et al.l d2005h and references 



therein), X-ray i Chernyakova et al.T 20061) and references therein), 
and TeV dAharonian et alj|2005h energy ranges. The orbital multi- 
wavelength variability pattern is determined by the details of the 
interaction of a relativistic pulsar wind with a strongly anisotropic 
wind of the companion Be star, composed of a fast, rarefied polar 
wind and a slow, dense equatorial decretion disk. The disk of the 
Be star in the PSR B 1259-63 system is believed to be tilted with 
respect to the orbital plane. While the inclination of the disk is not 
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constrained, the line of intersection of the disk plane and the orbital 
plane is known to be ori ented at about 90° with respect to th e major 
axis of the binary orbit dWex et al.lll998l : IWang et alj2004 and the 
pulsar passes through the disk twice per orbit. 

The unpulsed radio emission from the system appears approx- 
imately at the moment of the pre-periastron entrance of the pulsar 
into the equatorial Be star disk. Within several days unpulsed radio 
emission sharply rises to a peak, and then slightly decreases, as the 
pulsar passes through periastron. The second peak is reached dur- 
ing the second, post-perias tron disk crossing jjohnston et ail 19991 
l2005l : IConnors et al.l2002l) . 

ASCA observations of the PSR B 1 259-63 system in 1 994 and 
1995 have shown that, similar to the radio light cu rve, the X- 
ray light curve has two pe aks around the periastron l lKaspi et"aT1 
ll995MHiravama et all 1999b . The XMM-Newton observations of the 
source's orbital X-ray light curve, with a detailed monitoring of 
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Figure 1. Contour plot of the MOS2 XMM-Newton field of view for the X12 
observation in equatorial J2000 coordinates. A total of 4 contours were used 
with a square root scale between 5 counts per pixel (outer contour) and 500 
counts per pixel (inner contour). 

the period when pulsar approaches and enters the dense equatorial 
wind of the Be-star prior to periastron, have shown that the appear- 
ance of the unpulsed r adio emission is also accom panied by a sharp 
rise of the X-ray flux l lChernvakova et al. 112006b - The source spec- 
tral evolution has revealed an unexpected hardening of the source 
spectrum with the smallest value of the photon index T ~ 1.2, and 
a subsequent softening on the day scale as the pulsar moves deeper 
inside the disk. Unfortunately, in 2004 the source became invisible 
for XMM-Newton just after the entrance to the disk, so that the be- 
haviour of the source within the disk and during the second disk 
crossing remained unclear. In order to clarify this behaviour, we 
organised an intensive X-ray monitoring campaign during the 2007 
periastron passage. We have monitored the source with the Suzaku, 
XMM-Newton, Chandra and Swift satellites during a 2 months pe- 
riod, which covers both the pre- and post-periastron disk crossings. 
The results of this campaign are discussed below. 

This paper is organized as follows: in Section 2 we describe 
the details of the data analysis. The results are presented in Section 
3, and discussed in Section 4. 



2 OBSERVATIONS AND DATA ANALYSIS 

In 2007 we were able to organize an intensive X-ray monitor- 
ing campaign of the PSR B 1259-63 system with Suzaku, XMM- 
Newton, Chandra, and Swift. The list of observations is given in 
Table [T] where r denotes the number of days after the periastron 
passage (27 July 2007, MJD 54308) and is the true anomaly of 
the source. 

2.1 XMM-Newton observations 

In all XMM-Newton observations the source was observed with 
the European Photon Im aging Cameras (EPIC) M OS1, MOS2 
dPen Herder et al. 11200 lh and PN dStruder et ai]|200ll) detectors in 
the small window mode with a medium filter. The XMM-Newton 
Observation Data Files (ODFs) were obtained from the online Sci- 
ence Archiv^] and analyzed with the Science Analysis Software 

1 http://xmm.vilspa.esa.es/external/xmm_data_acc/xsa/index.shtml 



Table 1. Journal of 2007 of PSR B 1259-63. 



XMM-Newton observations 



Data 


Date 


MJD r <f> 


Exposure fjnosl 


f_mos2 


Set 




(days) (deg) 


(ks) 






Xll 
X12 
X13 


2007-07-08 
2007-07-16 
2007-08-17 


54289 -19 86.14 
54297 -11 112.26 
54329 +21 278.57 


9.34 
36.54 
6.35 


1.00!0;°> 


1 03+ 001 
1 05 +001 
1 05+ 001 


Chandra observations 


Chi 


2007-07-28 


54309 1 187.61 


4.68 






Ch2 


2007-08-06 


54318 10 248.87 


4.67 






Ch3 


2007-08-24 


54337 29 289.56 


3.15 






Ch4 


2007-09-18 


54362 54 308.83 


7.12 






Swift observations 


Swl 


2007-07-07 


54288.6 -20 84.1 


2.72 






Sw2 


2007-07-09 


54290.6 -18 88.3 


5.13 






Sw3 


2007-07-11 


54292.3 -16 92.6 


4.56 






Sw4 


2007-07-13 


54294.7 -14 99.5 


4.34 






Suzaku observations 


Data 


Date 


MJD r <f> 


Exposure 






Set 




(days) (deg) 


XlS(ks) 






Szl 


2007-07-07 


54288.6 -19.3 84.6 


21.9 






Sz2 


2007-07-09 


54290.7 -17.2 88.1 


19.5 






Sz3 


2007-07-11 


54292.6 -15.3 95.8 


22.7 






Sz4 


2007-07-13 


54294.7 -13.2 102.2 


22.9 






Sz5 


2007-07-23 


54304.3 -3.6 149.6 


19.7 






Sz6 


2007-08-03 


54315.3 7.4 230.0 


24.0 






Sz7 


2007-08-18 


54330.1 22.2 279.6 


20.5 






Sz8 


2007-09-05 


54348.2 40.3 300.4 


18.3 







(sas) v7.1.2. During the data cleaning, all time intervals in which 
the count rate in the energy band above 10 keV was higher than 1 
cts/sec for the PN detector and/or higher than 0.35 cts/sec for the 
MOS detectors, have been removed. In addition, we discarded the 
first kilosecond of X12 and the last two kiloseconds of X13 obser- 
vations, as these were affected by soft proton flares. 

The event lists for the spectral analysis were extracted from a 
15" radius circle at the source position for the Xll MOS1 obser- 
vation, and from a 22.5" radius circle for all other MOS and PN 
observations. We have performed spectral analysis by simultane- 
ously fitting the data of MOS l,MOS2 and PN instruments, leaving 
the inter-calibration factors between the instruments free. The val- 
ues of the MOS 1 and MOS2 inter-calibration factors relative to the 
PN are given in last two columns of Table[T] 

2.2 Chandra observations 

During our 2007 mo nitoring campaign, PS R B 1259-63 was ob- 
served with Chandra JWeisskopf et alJuOOOl) four times, see Table 
[T] The first two observations (Chi and Ch2) were done in the Faint 
mode with HETG grating. In order to reduce the pile-up, only 1/8 
of the S3 detector was in use, which reduces the nominal frame in- 
tegration time from 3.2 seconds to 0.4 seconds. The observations 
Ch3 and Ch4 were done in CC33_FAINT mode, because of the 
higher expected source flux level. This allowed us to obtain good 
quality results not affected by a pile-up. Similar to the first two 
observations, only the S3 detector was in use. In the data analysis 
we used the Chandra Interactive Analysis of Observations software 
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Figure 2. PSR B 1259-63 orbital evolution of (1-10 keV) light curve (top panel), spectral index (middle panel) and hydrogen column density (bottom panel), as 
seen with Suzaku,XMM-Newton,Swift and Chandra during the 2007 periastron passage along with the old XMM-Newton,BeppoSAX and ASCA observations. 
1-10 keV flux of the source is given in units of 10~" ergs cirT 2 s~' . Shaded area corresponds to the disk position proposed in Chemyakova et al. (2006). The 
solid line on the top panel is a fit with a Gaussian decay model (see text). To guide the eye we have connected old XMM-Newton data on the middle panel. 



Table 2. Spectral parameters of IGR J13020-6359 during 2007 Swift and 
XMM-Newton observations. 



Set 



Swl 
Xll 
Sw2 
Sw3 
Sw4 
X12 
X13 



10 



F(2-10 keV) 

-il -2 -! 
ergs cm _ s 


r 


10 22 cm" 2 


X 2 (dof) 


2 10+ 065 
-fl-3? 
2 21+" 1 ' 

z - -0.62 

1,90 95 

1 oo+?^f 
1.83 

t 4+0.06 

J - -0.09 

2 89 +01 ° 
z - -0.22 


1.00+ 118 
0.97°-^ 4 

, 9 ,+858 

1Z1 -0.I9 
1 44+0.27 

0.89+i 
0.87+0$ 


3.00 
2 27+0.I8 

2 56+^ : ^ 
ZJO -0.40 

3 j 4+0.69 

3 96 +( ™ 
233+8:^ 

9 O9+0.15 
Z - JZ -0.14 


21.969(11) 
171.08 (153) 
12.28 (22) 
7.28 (14) 
11.21 (11) 
534.78 (422) 
151.03 (150) 



* Due to the lack of statistics we fixed the value of N# in the Swl 
observation. 



package (CIAO ver. 3.4) and the CALDB version 3.4.0. The tool 
celldetect found no sources other than PSR B 1259-63 in the 
field of view. We have used the generated region file for the spec- 
tral extraction, and collected background events from a nearby re- 
gion of the same form and size. The resulting spectra were grouped 
to have at least 30 counts per energy bin. PSR B 1259-63 showed 
a moderate level of activity, with an average count rate of 0.21 
cts/s in Chi observations and 0.15 cts/s in Ch2 observation, which 
corresponds to a pile-up affection at the 1 5% and 10% levels, re- 
spectively0 It is worth to noteworthy that lKishishita et all d2009l) 
in their analysis of Chandra LS 5039 data found that even a rel- 



2 Pile-up affection level was calculated using PIMMS v3.9d tool, available 
at http://cxc.harvard.edu/toolkit/pimms.jsp Using this tool we assumed that 
emission can be described as an absorbed power law with parameters listed 
in TablefJ 
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Figure 3. Top paneUVSR B1259-63 orbital light curve (1-10 keV), as seen with SuzakuJCMM-Newton,Swift and Chandra during the 2007 periastron passage 
along with the old XMM-Newton,BeppoSAX and ASCA observations. 1-10 keV flux of the source is given in units of 10 _1 1 ergs cnT 2 s~' .middle panel: 2004 
HESS orbital light curve, bottom panel: Collection of historical 1.4 GHz radio light curve during different periastron passages. In order to compare the orbital 
behaviour at different wavelength we show with a black dash line the scaled 2004 radio curve on top and middle panels. Shaded area corresponds to the disk 
position proposed in Chernyakova et al. (2006). 



ative small pile-up fraction can result in somewhat harder photon 
indexes. Thus the absolute numbers found for Chi and Ch2 should 
be treated with some care. 



2.3 Swift observations 

The Swift/XKT i lGehrels et all |2004|) data were taken in photon 
mode with a 500 x 500 window size. We have processed all the data 
with standard procedures using the FTOOLS|jtask xrtpipeline 
(version 0.11.6 under the HEAsoft package 6.4.0). We have ex- 
tracted source events from a circular region with a radius of 20 
pixels (1 pixel a* 2.27"). To account for the background, we 
have extracted events from a nearby circle of the same radius. 



See http://heasarc.gsfc.nasa.gov/docs/software 



Due to the low count rate (less than 0.4 cts/s) no pile-up correc- 
tion was needed. The spectral data were rebinned with a mini- 
mum of 25 counts per energy bin for the x 1 fitting. The ancil- 
lary response file was generated with xrtmkarf, taking into ac- 
count vignetting, and the point-spread function corrections. In our 
analysis we have used the swxpcOtol2sQ_2Q010101vO10 . rmf 
spectral redistribution matrix for observations Swl, Sw2 and Sw3 
and the swxpc0tol2s6_2O010101vO10 . rmf spectral redistribu- 
tion matrix for the observation Sw4. 



2.4 Suzaku observations 

Suzaku has intensively monitored the PSR B 1259-63 2007 perias- 
tron passage. Observat ions were performed with the X-ray Imag- 
ing spectromet er (XIS:lKovama et al.l | |2007|) ) and the Hard X-ray 
detector (HXD: lTakahashi et al. 1 20071) ). These data were first pre- 
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sented in lUchivama etal] d2009h . For completeness we present the 
list of Suzaku observations in TableQ] 



3 RESULTS 

3.1 Imaging Analysis 

X-ray emission from PSR B 1259-63 has been clearly detected 
in all observations with all instruments. Apart from PSR B1259- 
63, another X-ray bin ary, IGR J 13020-6359 (iden tified with 2RXP 
J130159.6-635806 in IChernvakova et all l l2005h ). located 10 ar- 
cminutes away to the north-west, is detected in the MOS2 XMM- 
Newton and XRT Swift fields of view. Figure Q] shows the con- 
tour plot of XMM-Newton field of view for the X12 observation, 
in which both sources are visible. IGR J 13020-6359 was not de- 
tected in MOS1, because some of its CCDs were turned off during 
the observation. This source is found to be in a low activity state, 
with the flux comparable to the flux of PSR B 1259-63. Tabled 
summarizes the spectral parameters of IGR J 13020-6359 found in 
our observations. 



3.2 The X-ray light curve 

The upper panel of Fig. [2] shows the X-ray orbital light curve of 
PSR B 1259-63 system. The figure summarizes all the available 
X-ray data on the source in 1-10 keV energy range. The histori- 
cal data of XMM-Newton (XI - X10 ) and BeppoSAX points are 
taken fromlChernvakova et al. I d2006h , ASCA data are taken from 



Table 3. Spectral parameters for 2007 observations of PSR B 1259-63. 



iHiravama et al. 1 19991) . Observations made with different instru- 



ments at close orbital phases are consistent with each other, demon- 
strating satisfactory inter-calibration between different instruments: 
the points marking Szl and Sz2 and Xll and a Swift observation 
at the phase 9 ~ 85° are superimposed on each other. Some of our 
new observations are done at orbital phases close to the ones of his- 
torical ASCA observations (e.g. Chi). One can see that the system's 
orbital light curve does not exhibit strong orbit-to-orbit variations - 
all fluxes measured by ASCA are consistent with the ones measured 
13 years later. Our new data also confirm a local minimum at the 
phase 6 m 90°, first observed with XMM-Newton during the 2004 
periastron passage (X9 and X10 observations). 

Stability of the orbital light curve allows us to use old and 
new data simultaneously while analyzing the flux orbital evolution. 
The first (pre-periastron) entrance of the pulsar into the equatorial 
disk of the Be star is accompanied by a sharp rise of the X-ray flux 
(points X7 - X9), by a factor of six in seven days (Ad = 13°). This 
period of the sharp rise of the flux is followed by a period of vari- 
ability, most probably related to the interaction of the pulsar wind 
with the Be star disk. During this period, the X-ray flux decreases 
by a factor of 1.5 for two days (Ad = 7°) and then rises again for 
the subsequent two days (A6> = 6°). The decay of the flux after the 
exit of the pulsar from the disk is slower than the rise at the disk 
entrance. The flux decreases by a factor ~ 2.5 up to the moment of 
the post-periastron entrance to the disk. The second entrance to the 
disk is again accompanied by a sharp rise of the flux, by a factor 
of 3.5 in less than ten days (A9 < 30°). The second peak is fol- 
lowed by an adiabatic decay with the characteristic decay scale of 
30°, illustrated in Fig.[2]by solid curve, representing the best fit of 
the data with Gaussian decay model (f(6) = a * e ~ w ~ eo) ~ IA6 o + c, 
6 = 28O°,A0 O = 30°). 

For comparison we also show in Fig. [3] the TeV light curve 
of 2004 HESS observation dAharonian etai] |2005|) and radio 
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<P 


A 1 1 lu rvC V J 


r 


M 

n H 




Swl 
Szl 


84.1 
84.58 


2.58 + °36 
2.58±0.03 


1.64 ±0.02 


53+ - 07 

— U.Uo 

0.5 ± 0.02 


33 (35) 
375(354) 


Xll 

Sw2 
Sz2 


86.14 
88.3 
88.10 


2 - 76+ ool 

— U.U4 

9 4Q+0.25 
z - -0.22 
2.67±0.03 


| coq+0.010 
— U.UlU 

1 56 +0 - 07 
i.jo _ 0Q7 

1.58 ±0.03 


0.478 + 0-^ 
— y.uuo 

52+ - 05 
UJZ -0.05 

0.51 ± 0.03 


1490 (1575) 
63 (654) 
321(309) 


Sw3 
Sz3 


92.6 
95.80 


7 oq+0.27 

2.16±0.03 


i oc+0.09 
-0.08 

1.35 ±0.03 


n 40+0.06 
"•^-0.06 

0.51 ± 0.03 


41 (48) 
324(272) 


Sw4 
Sz4 


99.5 
102.17 


2 41 +0.27 

-0.24 

2.31 ±0.03 


i C7+0.08 

-0.08 

1.44 ±0.03 


56+ - 06 

-0.05 

0.52 ± 0.03 


42 (53) 
387(298) 


X12 
Sz5 


112.26 
149.61 


2 S4+ 02 
z ""-0.02 

1.9()±0.02 


1 530+ 0005 
1 """-0.005 

1.83 ±0.03 


(1 47X+0-003 

0.50 ± 0.03 


2792 (2543) 
212(249) 


Chi 


187.61 


1 24+ 006 
L - -0.07 


1 65 +0 - 07 


0.48 


16(30) 


Sz6 


230.02 


1.29±0.02 


1.73 ±0.03 


0.52 ± 0.03 


189(205) 


Ch2 


248.87 


0.93+0.06 


1,46+oos 


0.48 


23 (22) 


X13 
Sz7 


278.57 
279.61 


3 40+°' :o ^ 

J - -0.04 

3.44±0.04 


1 647+O-om 

" -0.009 

1.69 ±0.02 


0.480+^ 
0.50 ± 0.02 


1648 (1595) 
406(402) 


Ch3 

Sz8 


289.56 
300.35 


3 5Q+O.O6 
J - -0.07 

2.27±0.03 


1 47+0.03 
1.69 ±0.02 


0.48+™ 2 
0.46 ± 0.03 


187 (196) 
283(253) 


Ch4 


308.83 


1 43+0.03 

-0.03 


1 48+ 004 

1 -^°-0.04 




193 (200) 



In units of 10 11 erg cm 2 s 1 . 
* In units of 10 22 cm~ 2 . 



jjohnston et al.l 1 19991: IConnors et all 120021 : Ijohnston et all 120051) 
light curves of different periastron passages. It is seen that the gen- 
eral flux behaviour (rapid rise at the disk entrance, followed by a 
slow decay) in all energy bands is similar. To make this point more 
clear we add the scaled radio light curve of 2004 periastron passage 
to panels showing the X-ray and TeV data (dashed line). 



3.3 Spectral Analysis 

The spectral analysis was done with the NASA/GSFC XSPEC 
vl 1.3.2 software package. In Fig.|4]the folded and unfolded spectra 
of PSR B 1259-63 for X12, Sw3, Ch2, and Ch3 observations are 
shown. In order to make the figure clear, we have multiplied the 
observed spectra by the specified factors. 

A simple power law model with photoelectrical absorption de- 
scribes the data well in most observations, with no evidence for any 
line features. In Table[3]we present the results of the three parameter 
fits to the Suzaku, XMM-Newton, Chandra and Swift data in the 0.5 
- 10 keV energy range. The uncertainties are given at the lcr statis- 
tical level and do not include systematic uncertainties. The quality 
of the first two Chandra observations prevent the simultaneous de- 
termination of the spectral slope and absorption column density, so 
we decided to fix the latter to the value of Nh = 4.8 x 10 21 cm _1 , 
consistent with the value found in XMM-Newton observations. Or- 
bital evolution of the spectral parameters is shown in Figure [2] 
The high value of the reduced x 1 f° r Sz3 (1.19 for 272 d.o.f.) 
a nd Sz4 (1.30 for 298 d .o.f.) observations was studied in details 
in lUchivama etal] j2009h . It was found that for these observations 
broad band (0.6-50 keV) spectrum is much better fitted with a bro- 
ken power law model with a spectral break from Tj = 1.25 ± 0.04 
below Etr ~ 5 keV to T2 = 1.6 ± 0.05 above. 
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Figure 5. Light Curves of XI 1 (left), X12(middle) and xl3(right) observations in 0.3 - 2 keV (top panel), 2-10 keV (middle panel) along with hardness ratio 
(bottom panel). Time bin size is equal to Iks. Dashed line indicates the best fit of the data with a constant model. 



PSR B1259-63 




Energy, keV 

Figure 4. Folded(top panel) and unfolded (bottom panel) PSR B 1259-63 
spectra during the X12 (PN data), Sw3, Ch2 and Ch3 observations. In order 
to make the figure clearer we have multiplied the observed spectra by the 
specified factors. 



3.4 Timing Analysis 

Within one orbit PSR B 1259-63 is known to be variable at differ- 
ent time scales. Apart from the orbital time scale (several years), 
fast, day-scale, variability of the flux and spectral characteris- 
tic is observed during the periastron passage. Variability at time 
scales much shorter than the orbital time scale is observed in 
other TeV binaries. For example, the X-ray emission from the bi- 
nary LSI +61° 303 is characterized by a spectral energy distri- 
bution similar to the one of PSR B 1259-63 (which might point 
to the similar mechanism of the X- ray emission) wit h variability 
time scales of down to ~ 1 ksec dSidoli etalj|2006l) and possi- 
bly even shorter dSmithetal .1120081) . Such short time scale vari- 
ability c ould be related to the clumpy structure of the wind of the 
Be star dZdziarski et~aT]|2008l) . If so, then the clumpy structure of 
the Be star wind should also lead to the short-time scale variabil- 
ity of emission from the PSR B 1259-63 system. In order to study 



the variability on the short timescales, we have analyzed the light 
curves in the 0.3-2 keV and 2-10 kev energy ranges in individual 
X-ray observations. 

Figure [5] shows the soft and hard energy band light curves 
(Iks time bins) of the three 2008 XMM-Newton observations of the 
source. The lower panels of this figure show the hardness ratios. 

The source flux appears to be variable at the observation (~ 
10 ksec) time scale at least in the Xll and X12 observations. In 
the case of the observation Xll, the best fit of the soft (hard) band 
light curves with a constant flux, shown in the Figure|5] gives a^ 2 = 
55.88 Of 2 = 42.52) for 13 degrees of freedom. The probability that 
the soft (hard) band flux stayed constant over the entire observation 
is 3 x 10~ 7 (5 x 10~ 5 ), i.e. the constant flux hypothesis is ruled out 
at a ~ 6cr (~ 5cr) level. 

To study the details of variability at the ~ 10 ks time scale, 
we arranged a long ~ 50 ks continuous exposure during the ob- 
servation XI 2. The flux of the source was stable (1.4 cts/s in 2-10 
keV energy range) during the first ~ 30 ks of the observation and 
then during the following 10 ks grew up to a count rate of 1.6 cts/s. 
No sharp variations of the flux and hardness ratio on the short time 
scales have been observed. No significant variations of the flux are 
detected in the XMM-Newton observation X13: the scatter of the 
data points in the light curve shown in the right panel of Fig.[5]does 
not exceed the 3<x limit in the soft or the 2cr limit in the hard energy 
bands. 

A possible episode of faster (< 10 ks time scale) flux vari- 
ations was detected in the Xll observation, which is close to the 
moment of pulsar entrance into the disk and is in the middle of 
the period of strong spectral variability. In order to find the min- 
imal variability time scale in this observation we have applied the 
structure function analysis dSimonetti et aljl985h to the X-ray light 
curve. 

Fig. [6] shows the structure functions calculated for the soft 
band light curve of Xll observation. The error bars of the struc- 
ture function are estimated via the simulation of 10 4 light curves in 
which the values of the flux in each time bin are scattered around 
the measured values of the flux. The structure function deviates by 
3cr from the low plateau value at the time scale r - 3 ks, which 
means that the light curve is variable at this time scale. For com- 
parison we show, in the same figure, an example structure function 
of one of the simulated light curves, in which the flux variations are 
only due to the statistical uncertainty of the signal. One can see that 
the structure function of the simulated light curve is consistent with 
the lower plateau value in all the considered time range. 
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Figure 6. Red: structure function of the 0.3-2 keV band light curve of the 
observation XI 1. Dashed horizontal line shows the low plateau. Green: ex- 
ample of the structure function of a simulated light curve in which flux 
variations are due to the statistical scatter of the signal. 
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Figure 7. Light Curves of Sw3 observation in 0.3 - 2 keV (top panel), 2 - 
10 keV (middle panel) along with hardness ratio (bottom panel). Time bin 
size is equal to Iks. 



Evidence of the variability of the source on the several 
kiloseconds time scal e was also r eported in some of ROSAT 
dCominskv et al.ll 1994b and Suzaku dUchivama et alj|2009l) obser- 
vations. 

An example of the Swift light curve during the Sw3 observa- 
tion is shown in Figure UJ One can see that the source monitoring 
strategy was to take a set of ~ 1 ks long snap-shots separated by 
longer (several kiloseconds) time intervals. From Fig. UJ one can 
see that no significant variability at the time scale ~ 10 4 s is found 
in the Swift data. 



4 DISCUSSION 

The 2007 observational c ampaign confirms the previous finding 
dChernvakova et al. IbOOfj) that during the first, pre-periastron pas- 



sage of the disk, the X-ray spectrum of the source hardens on a 
day scale reaching a value of the photon index smaller than 1.5. 
During the 2004 periastron passage, a spectral state with photon 
index T p h - 1.2 was observed at the beginning of the rapid growth 
of the X-ray flux at the disk entrance (<p = 73°). Remarkable, the 
time of the subsequent spectral softening coincided with the time 
of a step-like increase of the hydrogen column density (0 = 83°), 
which can be related to the penetration of the pulsar in the dense 
equatorial disk. The observed hardening, along with the successive 
soften ing and sharp flux rise, was attributed bv lChernvakova et al. I 
fcOOrj) to the injection of high energy electrons at the disk entrance 
(e.g. due to the proton-proton collisions), or to a sharp decrease 
of the high energy electron's escape velocity accompanied with 
the modification of their spectrum by Coulomb losses. An alterna- 
tive explanation of the hard X-ray spectrum observed at the disk 
entrance is possible within a model in which X-ray emission is 
produced via the sync hrotron mechanism jKhangulvan et alj2007l : 
lUchivama e"tai]|2009l) . In this case hardening of the spectrum to 
photon indexes r p h < 1.5 on a day time scale can be achieved if the 
electron energy loss is dominated by inv erse Compton (IC) scatter- 
ing lo sses in the Klein-Nishina regime jKhangulvan & Aharonianl 
1 20051) . The Klein-Nishina regime of IC scattering becomes im- 
portant at energies above E e > m^c /(2.7kT,) ~ 30 GeV, where 
T, is the temperature of the Be star. Decrease of the IC cross- 
section at high energies leads to the decrease of the efficiency of 
the IC cooling of electrons and, as a result, to the hardening of the 
electron spectrum. Electrons with a power law energy distribution 
N e oc E~ r cooled in the Klein-Nishina regime form a spectrum 

N e (E e ) oc E; r+[ X [in ^ - - 0. 6472]'', much hard er than 

N e (E e ) oc £"~ <r+l> in Thompson limit dBlumenthal&Gouldl 1 197ph . 
Thus the resulting synchrotron spectrum is proportional to e~ r ' 2 
and can be harder than 1.5. 

The 2007 data did not cover the period of the entrance to the 
disk, so it is not possible to check if the spectral hardening proceed- 
ing the flux growth repeats from orbit to orbit. Instead, the photon 
index F 



P h 



1.3 was observed in 2007 at almost the maximum of 
the X-ray flare associated to the pre-periastron disk passage. More 
precisely, the time of the hardening and subsequent softening of the 
spectrum in the 2007 data coincides with a local "dip" in the X-ray 
light curve during the broad flare associated to the disk passage. 
No strong variations of the hydrogen column density are notice- 
able at the moment of the hardening of the spectrum. Detection of 
the spectral break in the X-ray spectrum of the source at E ai 5 keV 



in Suzaku observation, which reveals the hard state with r. 



pi, 



1.3 



dUchivama et alj|2009h . provides an additional clue for the under- 
standing of the nature of the observed spectral hardening. 

These two models of spectral hardening could be readily dis- 
tinguished via multi-wavelength observations. In the case of the IC 
mechanism of the X-ray emission, the energies of electrons respon- 
sible for the X-ray emission are in the range of ~ 10 MeV, while in 
the synchrotron model the X-ray emitting electrons have multi-TeV 
energies. A necessary condition for the hardening of the spectrum 
beyond T = 1.5 in the synchrotron model is that the energy losses 
of multi-TeV electrons are dominated by the IC loss. This means 
that the expected TeV-band IC luminosity of the source at the mo- 
ment of the spectral hardening is larger than the X-ray luminosity. 
Unfortunately, the absence of the TeV observations simultaneous 
with the X-ray observations in 2007 does not allow us to distin- 
guish between the two models. 

If the hard X-ray spectrum can be explained by Coulomb 
losses of the ~ 10 MeV electrons, then the break at the 
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Figure 8. Spectral energy distribution of PSR B 1259-63 in the model 
with the power law injection spectrum of electrons with spectral index 
r e = 2. 2. Top panel shows the initial electron injection spectrum with the 
dashed line and the spectrum formed in result of cooling and escape of 
electrons to the distance D ~ 10 14 cm (solid line). Bottom panel shows the 
spectra of synchrotron (thin solid line), inverse Compton (dashed thin line) 
and bremsstrahlung emission (dotted thin line) as well as the overall spec- 
trum modified by the pair production in the photon field of the Be star (thick 
solid line). Blue and red lines at the radio part of the spectrum illustrate the 
slope and possible flux range of the observed radio emission 



~ 5 keV energy, observed by Suzaku, has to be ascribed to 
the so-called "Coulomb" break in the electron spectrum. Such 
a break appears at the energy E c , at which the rate of the 
Coulomb energy loss is equal to the IC energy loss rate: E c ^ 

15 [D/10 13 cm] [;! t ,/10 8 cirr 3 ]' /2 MeV, where n e is the stellar wind 
density at a distance D (10 13 cm is a characteristic separation 
of the companions at orbital phases close to periastron). Elec- 
trons with the energy E c produce IC emission in the energy band 
e c ~ 4[T,/3 x 10 4 K] [£ c /10MeV] 2 keV. Measurement of the 
break in the X-ray spectrum of the source enables an estimate of 
the density of the medium 



n e (D =a 10 13 cm) = 10 s cm" 3 



(1) 



This estimate is consistent with the one expected for an equatorial 
disk with the radial density profile n e (D) = noC-D/S*) -3 ' 5 if the disk 
density close to the surface of the star is no ~ 10 12 cirr 3 . 

If the X-ray emission is produced via the IC mechanism, the 
break in the X-ray spectrum can be also related to the existence of a 
low-energy cut-off in the electron spectrum. Indeed, the IC cooling 
time of the X-ray emitting electrons of the energy E is 



tlC(T) 



3mn 2 c 4 D 2 



VtL+E 
6x 10 5 



10 38 erg s" 1 



D 



10 13 cm 



lOMeV 



(2) 



where L+ is the luminosity of the Be star, m e is the elec- 
tron mass and a T is the Thomson cross-section. The IC cool- 



ing time is comparable to the escape time f csc a D/V — 

10 6 [o/10 13 cm] [v/10 7 cm/s] s if the escape velocity is V ~ 

10 7 cm/s. If the IC emitting electrons are initially injected at ener- 
gies much la rger than 10 MeV (e.g. as a resul t of the proton-proton 
interactions, iNeronov & Chernvakoval d2007t) ). they would not be 
able to cool to energies below ~ 10 MeV, which can explain the 
deficiency of the IC emission at the energies below the ~ 5 keV 
break energy. 

Finally, the shape of electron spectrum below 10 MeV can be 
also affected by adiabatic cooling. Equating the adiabatic cooling 
rate dEJdt = (EJR)dR/dt = nVE e /(2D% where R ~ D"> 2 is the 
radius of the synchrotron/IC emitting bubble and V is the escape 
velocity, to the IC cooling rate, dE e /dt = (4/3)o" r c[/ p h(£ e /m e c 2 ) 2 , 
where U P h = L*/4nD 2 c is the radiation energy density, one finds 
the adiabatic break energy 



10 MeV 



D 



10 13 cm 



V 



10 38 erg/s 



eV (3) 



10 7 cm/s 

If the emission is produced at the distances D ~ 10 13 cm, the ob- 
servation of the cooling break at the energy e =i 5 keV is consistent 
with the assumption that the break is produced by the influence of 
the adiabatic losses if the escape velocity is V ~ 10 7 cm/s, which is 
in agreement with the estimates of the velocity of the stellar wind in 
the considered range of distances ( I Waters. Cote & Lamersfl987T) . 

The results of numerical modeling of the broad band spec- 
trum of the source within the IC scenario of the X-ray emission are 
shown in Fig. [8] We assume that the high energy electrons are in- 
jected in the synchrotron/IC emitting bubble, which escapes with 
a speed equal to the speed of the stellar wind. The spectrum of 
the high-energy electrons, shown by the thick solid line in the up- 
per panel of Fig. [8] is formed as a result of cooling of electrons 
due to the IC, synchrotron, bremsstrahlung and Coulomb energy 
losses during their escape from the system. The injection of elec- 
trons is assumed to happen at a distance D = 3 x 10 12 cm. The 
injection spectrum of the electrons, shown by the dashed line in 
the upper panel of the fig ure is assumed t o be a power law with 
spectral index r in j = 2.2 dKirk et aL I l2000h . The magnetic field is 
assumed to be equal to So = 0. 1 G at the initial injection distance. 
In the calculations presented in Fig. [8] we have used the angle- 
averaged IC cross-section to calculate the IC energy loss and emis- 
sion spectrum. The account of the anisotropy of the IC emission 
is i mportant for the accurate ca l culati o ns of the source light curv e 
(see IChernvakova & Illarionovl dl999h : iKhangulvan et all fe008h ). 
but requires a detailed 3D model and is not addressed here for sim- 
plicity. The cooled spectrum of electrons is shown on the top panel 
of Fig. [8] with a solid line. Coulomb losses affect the electron spec- 
trum at the lower end, leading to the formation of the first bump. 
The second bump arises due to the decrease of the cross-section 
after the transition from Thomson to the Klein-Nishina regime, 
which results in spectral hardening above 10 u eV, and softening 
above 10 12 eV because of the dominance of the synchrotron losses 
at those energies. 

An alternative model, in which the observed X-rays are pro- 
duced via the synchrotron mechanism is illustrated in Fig. [9] For 
the calculation we have used the same numerical code as for the 
calculation shown in Fig. [8] The main differences between the two 
cases are different injection spectra and the assumed initial val- 
ues of the magnetic field. The magnetic field at the initial distance 
D = 3 x 10 12 cm is assumed to be Bo = 1 G. The electron injection 
spectrum, shown by the dashed line in the upper panel of Fig.|9]is 
assumed to be sharply peaked at energies about 10 13 eV. If the typi- 
cal initial electron injection energies are above 10 TeV, the shape of 
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Figure 9. Spectral energy distribution of PSR B 1259-63 in the model with 
electron injection at very high energies. Notations are the same as in Fig. [8] 



the cooled electron spectrum at the energies below 10 TeV is com- 
pletely determined by the cooling effects and is almost independent 
of the details of the initial electron injection spectrum. 

The low-energy cut-off in the cooled electron spectrum at 
E ~ 10 8 eV arises because at energies below 10 s eV the radia- 
tive cooling time of electrons becomes longer than the escape time 
from the system. The spectrum of electrons is hard in the 10 GeV - 
1 TeV energy interval, in which the electron energy losses are dom- 
inated by IC scattering proceeding in the Klein-Nishina regime. A 
sharp break in the spectrum at TeV energies is related to the fact 
that the cooling rate of the multi-TeV electrons is determined by 
the synchrotron, rather than IC energy losses. The break energy 
can be determined by equating the synchrotron loss time 



6nm 2 e c 3 
t s = — - 4 x 10 

cr T B 2 E, 

to the IC loss time in the Klein-Nishina regime, 



1G 


2 


1 TeV 


~B 







(4) 



\6E e D 2 h 



o- T (m e ckT i ,R < ,)' 1 
s\ E e If D 



In" 



0.55E e kT* 



10 



lTeVjll0 13 cm 



This gives 



0.5 



B 



D 



-1/2 



TeV 



(5) 



(6) 



.1 GJLlO 13 cm 

The synchrotron emission produced by electrons with the energy 
£t, r is emitted at the energy 

etiBE: 



10 



B 



1 GJL0.5TeV 



keV 



(7) 



which is close to the observed break energy e 5 keV in the Suzaku 
spectrum. 

The two different mechanisms of the X-ray emission (IC or 
synchrotron) could also be distinguished by the difference in the 



variability properties of the X-ray signal expected in the two mod- 
els. In principle, much faster variability is expected in the syn- 
chrotron model of X-ray emission. Within this model, the typical 
variability time scale is set up by the synchrotron and IC cooling 
times of the multi-TeV electrons. This time scale could be as short 
as < 1 ks, if the magnetic field in the emission region is strong 
enough. On the other hand, within the IC model of X-ray emission, 
the characteristic variability time scales are set up by the rate of the 
Coulomb and/or IC and/or adiabatic energy losses of ~ 10 MeV 
electrons, which are normally much longer than ~ 1 ks. The detec- 
tion of the variability in the XI 1 observation of XMM-Newton at 
the time scale of ~ 3 ks thus provides an argument in favor of the 
synchrotron model of the X-ray emission. However, fast variability 
in the XI 1 observation is detected at a 3<x level with the help of 
the structure function analysis. Further observations are necessary 
to firmly establish the presence/absence of the fast variability of the 
X-ray emission. 

To summarize, we find that the X-ray data alone do not al- 
low us to distinguish between the synchrotron and IC origin of the 
X-ray emission from the source. The observation of the harden- 
ing of the spectrum below T = 1.5 during the pre-periastron disk 
passage gives an important clue about the X-ray emission mecha- 
nism. However, the origin of the observed spectral hardening can 
be clarified with the help of the simultaneous TeV observations. If 
the observed X rays have an IC origin, then the observed hardening 
during the drop of the flux is primarly connected to the hardening 
of electron spectrum below ~ 10 MeV (due to the Coulomb losses, 
or because of the escape of electrons from the emission region), so 
that no tight correlation between the X-ray spectral evolution and 
the TeV energy band emission is expected. On the other hand, in 
the case of synchrotron origin of the observed X-rays, the spectral 
hardening can be produced if the electron cooling is dominated by 
the IC energy loss in the Klein-Nishina regime. This implies that 
the IC flux from the system in the very-high-energy band at the 
moment of the spectral hardening should dominate over the X-ray 
flux. GeV data can give us another possibility to distinguish be- 
tween the models, as in this region the predicted flux value is very 
different in IC and synchrotron model. Hopefully Fermi Gamma- 
ray telescope observations of the next PSR B 1259-63 periastron 
passage will give us a clue to the distinguish between the models. 
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